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Heat capacity and susceptibility measurements have been performed on the diluted spin dimer 
compound Ba3(Mni_i:Va;)208. The parent compound Ba3Mn208 is a spin dimer system based on 
pairs of antiferromagnetically coupled S — 1, 3d^ Mn^^ ions such that the zero field groundstate 
is a product of singlets. Substitution of non- magnetic S ^ 0, 3d° ions leads to an interacting 
network of unpaired Mn moments, the low temperature properties of which are explored in the 
limit of small concentrations, < x < 0.05. The zero-field heat capacity of this diluted system 
reveals a progressive removal of magnetic entropy over an extended range of temperatures, with no 
evidence for a phase transition. The concentration dependence does not conform to expectations 
for a spin glass state. Rather, the data suggest a low temperature random singlet phase, reflecting 
the hierarchy of exchange energies found in this system. 

PACS numbers: 75.10.Jm, 75.50.Lk, 75.30.Hx, 75.40.-s 



I. INTRODUCTION 

Randomness can lead to intriguing magnetic states not 
typically available to perfectly ordered systems. The 
archetypal example is the spin glass state, found for a 
wide variety of disordered materials with either site or 
bond randomness^. In contrast, for gapped systems with 
a singlet ground state, substitution of non-magnetic ele- 
ments can introduce local moments, ultimately leading to 
long range magnetic order due to the effective interaction 
mediated by the background singlet state - one manifes- 
tation of "order by disorder" (OBD). For example, both 
the Spin-Peierls compound CuGcOs and the spin ladder 
compound SrCu203 have gapped ground states without 
long range order; however diluting either system with a 
small amount of non-magnetic Zn or Mg ions onto the 
S ~ ^ Cn site induces antiferromagnetic order—"—. It is 
not clear whether this effect should be uniformly antici- 
pated for all gapped systems, motivating both theoretical 
and experimental interest aimed at exploring the prop- 
erties of disordered spin-gap materials. 

The broad category of spin dimer compounds provides 
a simple means to study the effect of non-magnetic sub- 
stitution on a singlet ground state. Comprising pairs of 
spins with a dominant antiferromagnetic nearest neigh- 
bor exchange, the ground state is a product of singlets. 
The gap to excited triplet states can be closed by an 
applied field, leading to canted XY antiferromagnetic 
order—. Recently, quantum Monte Carlo (QMC) sim- 
ulations have been performed for the specific case of 
spin dimers arranged on a square lattice with antifer- 
romagnetic nearest neighbor and next nearest neighbor 
exchange^ii. These calculations revealed that for an ap- 
propriate range of concentrations, substitution of non- 
magnetic impurities leads to long range order in zero 
magnetic field. The predicted wave-vector is the same 
as that found for the stoichiometric parent compound 
subjected to fields above the critical field. It would be 



highly desirable to experimentally test whether such an 
OBD state is found for a real material conforming to the 
simple effective spin Hamiltonian used in this calculation. 
Unfortunately, there are not currently any suitable can- 
didate materials that match these requirements for which 
substitution of non-magnetic species is possible over an 
appreciable range of concentrations. Conversely, it would 
be equally interesting to see how geometric frustration af- 
fects the stability of the OBD phase. Here, unfortunately, 
quantum Monte Carlo simulations are prohibitively dif- 
ficult due to the frustration induced sign problem, and 
we must resort to experiment to provide insight. In this 
case, though, we are much more fortunate in that there 
are several candidate materials to which we can turn. 
One such material is Ba3(Mni_2:Va;)208. In this pa- 
per, we present the results of an initial survey of the low 
temperature properties of this material via heat capac- 
ity measurements. We find no evidence for a sharp phase 
transition into an ordered state down to our base temper- 
ature of 50 mK. Rather, the magnetic entropy is found 
to be removed over an extended range of temperatures, 
nominally independent of impurity concentration for the 
range of compositions studied. These results do not con- 
form to expectations for a canonical spin glass, leading to 
the possibility that Ba3(Mni_a;Va;)208 manifests instead 
a random singlet phase at low temperature. 

Ba3Mn208 is a spin dimer compound composed of 
vertical pairs of Mn^"*" ions arranged on triangular lay- 
ers (Fig. [T^)}^"i£. Antiferromagnetic exchange between 
the S = 1, 3(P Mn^+ ions forming the dimer leads to 
a zero field ground state which is a product of singlets 
with excited triplet and quintuplet states. Interactions 
between dimers broadens the excited states, and appli- 
cation of a magnetic field can close the gap to the ex- 
cited states, leading to three distinct ordered states at 
high fieldsii"— . The interactions of this system have 
been determined through a combination of inelastic neu- 
tron scattering (INS), electron paramagnetic resonance 
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FIG. 1: (Color online) (a) Schematic diagram showing 
the transition metal sublattice of the dimer compound 
Ba3(Mni„3;Va;)20g. Mn ions are shown in red. Dominant 
antiferromagnetic exchange Jo between pairs of Mn moments 
in the undiluted compound (x—O) leads to a singlet ground 
state. In this diagram a single Mn ion has been substituted 
by a non- magnetic V ion (green). The Mn ion vertically be- 
low the V ion is now no longer paired and hence contributes 
a local magnetic moment. Distances between the V ion and 
its nearest neighbors ro - are labeled, (b) Probability, as- 
suming no clustering, of finding at least one V ion within a 
distance r from the substituted ion, shown for the specific 
vanadium concentrations (x) studied, and for the specific dis- 
tances ro — T/f^ (c) Probability of finding the nearest neigh- 
boring V ion a distance r from the substituted ion^i. 



(EPR) and thermodynamic measurements of the phase 
diagram. INS studies of Ba3Mn208 revealed a spin gap 
of A = 1.05 meV and a dominant exchange within a 
dimer (along the ro direction) as Jq = 1.61 meV— li^. 
These studies also determined the nearest and next near- 
est out-of-plane interdimer exchanges, Ji = 0.118(2) meV 
and Ji = 0.037(2) mcV; and the dominant in-plane in- 
terdimer exchanges J2 ~ J3 = 0.1136(7) meV (ri - r4 
in Fig. [T])(a). Measurements of the critical fields then 
yield estimates for J2 and J3 of 0.256 meV and 0.142 
meV, respectively^. Finally, EPR measurements of the 
Mn^+ moments in Ba3(Mni_a;Va;)208 for nominal com- 
position of a; = 0.75 revealed a nearly isotropic g-tensor, 
with gaa = 1-96 and gcc ~ 1.97, and an easy axis sin- 
gle ion anisotropy D = -0.024 meV— . Similar measure- 
ments in the pure Ba3Mn208 compound revealed a zero 
field splitting of the triplet states characterized by D =- 
0.032 meVi^. 

Ba3V208 is isostructural to Ba3Mn208. However in 
contrast to the strongly magnetic Mn^+ ion, V'"'"'' cor- 
responds to a 3^*^ electron configuration and hence is 
non-magnetic. Partial substitution of V in Ba3Mn2 08 
therefore leads to unpaired Mn moments (Fig. [Ua))^^. 
For the relevant V concentrations studied here, the data 
do not suggest clustering of impurities, in which case, for 
the highest dilution studied, x = 0.046, the probability 
is 0.64 to find at least one neighboring V ion within 8 A 



of any given V impurity (Fig. [Ifb)). Furthermore, the 
most likely distance between nearest neighboring V ions 
is the r2 pairing (Fig. [H^c)). 



II. EXPERIMENTAL METHODS 

Single crystals of Ba3(Mni„xVx)208 (0 < a; < 1) were 
grown from a NaOH flux, following the procedure that 
we have previously described for Ba3Mn20aii. Poly- 
crystalline precursor was made by mixing BaO, Mn203 
and V2O5 in appropriate ratios. Electron microprobe 
measurements using Ba3Mn2 08 and Ba3V2 08 standards 
were used to determine the vanadium content of the sin- 
gle crystals. Measured V concentrations were close to 
nominal values (Fig. [2jc)). Data are presented here 
for samples with V concentrations 0.009(1), 0.020(1), 
0.034(1), 0.046(1) and 0.980(1). Uncertainties reflect 
the standard deviation between multiple measurements 
performed at different locations for individual crystals. 
Systematic uncertainties for these low concentrations are 
likely slightly larger. Crystals of Ba3Mn208 have an in- 
tense dark green color. Crystals of Ba3V2 08 are perfectly 
clear due to the absence of 3d electrons. All samples were 
found to be insulating regardless of vanadium concentra- 
tion. 

Low-field susceptibility measurements were performed 
using a commercial Quantum Design MPMS XL SQUID 
magnetometer for fields of 1000 Oe applied perpendicular 
to the c-axis. Heat capacity (Cp) studies were performed 
with a Quantum Design physical properties measure- 
ment system (PPMS) using standard thermal relaxation- 
time calorimetry. These measurements were performed 
in temperatures down to 50 mK and fields up to 0.5 T 
parallel to the c axis. 



III. RESULTS 

The low field susceptibility of Ba3(Mni_2;V2;)208 for 
< a; < 0.046 is shown in Fig. [Ifa). Results were fit 
to a dimer model including a mean field correction, a 
term corresponding to Curie- Weiss paramagnetic behav- 
ior and a temperature independent term, according to 
the previously described model^i^i: 
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Here A = 3 [Ji + J4 -I- 2 ( J2 -f J3)] /NAg'^fJ^^ is a mean 
field correction to account for interdimer exchange. Na 
is Avogadro's number, /3 = l/fc^T and a is the number 
of dimers per mole (where 1 mole refers to the formula 
unit Ba3(Mni_a;V2;)208, such that a = 1 for a; = 0). 
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FIG. 2: (Color online) (a) Susceptibility as a function of 
temperature for single crystals of Ba3(Mni_a;Va;)208 (com- 
positions listed in legend). Successive curves are offset by 
0.004 emu/mol for clarity. Data are fit (red lines) by a dimer 
model, described in the main text, including a mean- field cor- 
rection to account for interdimer interactions, a temperature- 
independent term, and a Curie- Weiss term to account for the 
unpaired magnetic impurities introduced by V-substitution. 
Values for the exchange constants were held fixed at the 
values obtained for the stoichiometric compound (black cir- 
cles), (b) Curie constant, C (left axis), and the associ- 
ated effective moment per V, /^e/// V (right axis). Note 
that V^"'' is non-magnetic, and that the moment arises from 
the unpaired Mn^^ 5'=1 spin on the broken dimers. Dot- 
ted horizon tal line s hows the anticipated effective moment 
Me// = g^/S{S~vV)^^B = 2.83 /is for Mn'^+ (c) Measured 
versus nominal composition of the single crystals used in this 
study. 



Recent INS studies found a negligible concentration de- 
pendence of both the spin gap A and the triplet band- 
width for X < 0.05^. Therefore values of Jo = 16.42 
K, Ji + J4 + 2 ( J2 + J3) = 5.31 K and g = 2.07 de- 
termined from the fit of the undiluted Ba3Mn2 08 mea- 
surement (solid black circles) were held fixed for the 
fits of susceptibility data for the other samples. These 
fixed values are in approximate agreement with the val- 
ues measured from INS and EPR, of Jq = 18.78 K, 
Ji + J4 + 2 ( J2 + J3) = 11-03 K and g = 1.96, up to 
the inherent limitations of this fit. The Curie constant 
C, extracted from these fits, is plotted against the left 
axis in Fig. [2{b), and varies essentially linearly with V 
concentration. The effective moment per mole of V, plot- 



ted against the right axis of Fig. [HJb), is therefore inde- 
pendent of V concentration within the uncertainty of the 
measurement. Observed values are close to those antic- 
ipated for Mn^''^, implying that each V impurity results 
in a single unpaired Mn^+ spin. That is, substitution of 

V docs not appear to result in clustering, but rather, as 
anticipated in Fig. [TJc), the number of doubly-broken 
dimers is very small, and the majority of V impurities 
occupy half of a dimer site paired with a magnetic Mn^"*" 
ion. The fit parameter a nominally scales with (l-x). 
However, because the change in a is small relative to its 
absolute value for the low concentrations studied here, 
this provides a much less sensitive measure of the num- 
ber of unbroken dimers than does the Curie susceptibility. 
The temperature independent background term xo was 
motivated by Langevin diamagnetism, but fit parameters 
were found to be weakly paramagnetic (xo — 3±2 x 10~^ 
emu/mol for the parent compound), most likely refiecting 
inaccuracies associated with the dimer model, and specif- 
ically the relatively crude mean-field approximation used 
to account for interdimer exchange interactions. Uncer- 
tainty in fit values of xo were larger than any variation 
with composition. 

Heat capacity measurements were performed down to 
50 mK using a dilution refrigerator. Data are shown 
in Fig. [ll^a) for the temperature interval 50 to 500 mK. 
The stoichiometric parent compound, Ba3Mn208 (x = 0, 
black filled circles) exhibits only a very small heat capac- 
ity in this temperature interval due to the much larger 
energy scale of the spin gap. A small temperature de- 
pendence due to the phonon contribution is barely dis- 
cernable, but can be seen more clearly when similar data 
for Ba3V2 08 are plotted over a wider temperature inter- 
val in Fig. [Sl^b). In contrast, samples with x > re- 
veal dramatically different behavior. The heat capacity 
is approximately two orders of magnitude larger, and ex- 
hibits a broad maximum centered at approximately 110 
mK. For a;=0.009 and 0.020 there is a slight indication of 
an additional upturn in the heat capacity at the lowest 
temperatures. This is not observed for the higher two 

V concentrations, a;=0.034 and 0.046, which also exhibit 
a progressive broadening of the feature at 110 mK, and, 
for a;=0.046, a gradual increase in the heat capacity from 
200 mK to 500 mK. Significantly, none of the concentra- 
tions studied exhibit a sharp anomaly characteristic of a 
phase transition. 

The same data shown in panel (a) of Fig. [3] are replot- 
ted as Cp/T normalized by the vanadium concentration 
2a;, in panel (b). Scaled in this way, it is clear that Cp 
approximately scales with x. The change in entropy, /S.S 
over this temperature window can be estimated from the 
area under these curves, and is plotted against the right 
axis of panel (b), also normalized per mole of V. The 
small phonon contribution has negligible effect on this 
estimate over this temperature range, as can be readily 
appreciated by inspection of the data in panel (a). In- 
spection of panels (a) and (b) of Fig. [3]reveals a progres- 
sive change in the T-depcndence of the heat capacity and 
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FIG. 4; (Color online) Heat capacity of 

Baa (Mno. 980 Vo. 020)208 in magnetic fields of T, 0.25 
T and 0.5 T, shown as red, blue and green points, respec- 
tively. Lines show the theoretical heat capacity for the 
Schottky anomaly associated with an isolated S — 1 spin 
with single ion anisotropy D = -0.024 meV in the same 
fields. Inset shows the appropriate energy spectrum for this 
calculation. Vertical lines indicate the three different fields 
for which measurements were performed. 



FIG. 3: (Color online) (a) Low temperature heat capacity of 
Ba3(Mni„3;Va;)208 in zero field. Substitution of vanadium 
leads to a substantial increase in the total heat capacity, (b) 
Left axis: The same data, plotted as Cp/T and scaled by V 
concentration, revealing the evolution in the functional form 
of the heat capacity. Right Axis: Change in entropy AS 
between 50 mK and 500 mK associated with the total heat 
capacity. Despite subtle changes in the functional form of 
the heat capacity as the vanadium concentration is increased, 
the integrated entropy per mole of V over this temperature 
window remains remarkably similar. 



entropy as the vanadium concentration increases. For the 
higher vanadium concentrations, relatively more entropy 
is removed per increment of temperature at high tem- 
perature than for lower concentrations, and vice versa at 
the lowest temperatures. This concentration dependence 
is, however, rather subtle, and the most remarkable as- 
pect of the data is that the total entropy removed over 
this interval is almost identical for all four compositions, 
despite a factor of 5 difference in the vanadium concen- 
tration. This value [~ 6 J/K mole (V)] is approximately 
60% of the total magnetic entropy Stotai = i?ln(3) (where 
R is the molar gas constant) associated with the unpaired 
S=l magnetic Mn ions induced by the vanadium substi- 
tution. 



The heat capacity data for the V-substituted samples 
shown in Fig. |3l[a) reveal a broad maximum centered 
at approximately 110 mK for all four compositions. In 
the absence of any other interactions, uniaxial single ion 
anisotropy, represented by a term D {S^)^ in the spin 
Hamiltonian, will split the ~ Q and S"^ = ±1 triplet 
states of unpaired Mn moments. A calculation of the 
expected contribution to the specific heat in zero field 
is shown in Fig. 2] for the specific case oi x ~ 0.020 
(red line), using the value oi D = -0.024 meV previ- 
ously determined from EPR measurements of diluted 
Ba3(Mni_a;Vi:)208. The maximum value occurs at a 
temperature very close to that observed in experiment 
(solid red circles), indicating that the experimentally ob- 
served features are likely related to Schottky behavior. 
To further test this hypothesis, the heat capacity was also 
measured in applied fields of 0.25 T and 0.5 T (green and 
blue data points in Fig. |4|) and compared with calculated 
values (blue and green lines) . The applied field splits the 
S''- — ±.\ states (inset to Fig. S]), resulting in broaden- 
ing of the Schottky anomaly and a shift in the maximum 
value to higher temperatures. The experimental data re- 
veal very similar behavior, confirming that this feature is 
closely related to the Schottky behavior anticipated for 
unpaired Mn spins. 

The correspondence between the calculated Schottky 
anomaly and the experimentally determined heat capac- 
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ity is, however, not perfect. Inspection of Fig. ^reveals a 
significant discrepancy between the magnitude of the the- 
oretical curves and the experimental data, even though 
the temperature at which the maxima occur agrees. This 
difference is particularly striking for the T data, for 
which the measured heat capacity is uniformly larger 
than the theoretical prediction, indicating that the mag- 
netic entropy associated with the full triplet S = 1 state 
is progressively removed over a fairly broad tempera- 
ture interval. For non-interacting moments, the doublet 
ground state should remain unsplit. Hence the antici- 
pated entropy associated with the Schottky anomaly in 
zero field is S'^ = 2a;[i?ln(3) - i?ln(2)] = 2a:i?ln(3/2) 
(where 2x is the number of moles of free Mn ions). In 
contrast, the measured magnetic entropy between 50 mK 
and 500 mK (Fig. |3Jb)) is considerably larger than this 
value. Comparison of the measured heat capacity in 
zero field with the calculated Schottky behavior indicates 
that much of this difference occurs at a temperature con- 
siderably above the energy scale set by the single ion 
anisotropy. Higher order terms in the crystal field ex- 
pansion would only lead to splitting on a lower energy 
scale than the leading axial term and cannot account for 
this difference. Nor does it seem likely that the progres- 
sive removal of entropy above 0.2 K is due to a spread 
of values of D since the total entropy would still rise to 
2.Ti?ln(3/2), and also since the anomaly centered at 110 
mK is not especially broadened. The data therefore in- 
dicate that interactions between the free moments play 
a significant role. Indeed, the magnitude of the Schot- 
tky anomaly itself appears to be smaller than anticipated 
for the given concentration, but superimposed on top of 
a large and only weakly temperature dependent back- 
ground, which must then arise from these interactions. 
This is also borne out by measurements made in applied 
field. 

In the presence of a magnetic field the 5^ = ±1 doublet 
is split, leading to an overall increase in the magnitude 
of the calculated heat capacity (blue and green lines in 
Fig. S]). For 0.25 T the measured heat capacity exceeds 
the calculated value for temperatures greater than 0.4 
K, indicating the removal of magnetic entropy at higher 
temperatures, similar to the zero field data. The area 
under the calculated blue and green curves is the full 
magnetic entropy of the free S = 1 spins, corresponding 
to = 2a::i?ln(3). If some of this magnetic entropy is re- 
moved at higher temperature due to interactions between 
Mn moments, then there will be less magnetic entropy 
available at low temperatures, which is presumably why 
the measured heat capacity in 0.25 T falls below the cal- 
culated curve for temperatures below 0.4 K. Data for 0.5 
T appear to follow the same general form: more magnetic 
entropy is removed at higher temperature than would be 
anticipated for isolated Mn ions, reducing the magnitude 
of the low temperature Schottky anomaly. 

Measurements of the heat capacity for x — 0.020 were 
extended to higher temperatures in order to integrate 
the magnetic entropy over a wider temperature window 
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FIG. 5: (Color online) a) Heat capacity of 
Ba3(Mno. 980 Vo. 020)208, Ba3Mn208 and the difference 
at T shown in black circles, red squares and the blue line, 
respectively. The change in entropy calculated from the 
difference is plotted as a black dotted line along the right 
axis, b) Heat capacity of Ba3(Mno.o2oVo.9so)208, Ba3V208 
and the difference at T shown in black circles, red squares 
and the blue line, respectively. The change in entropy 
calculated from the difference is plotted as a black dotted 
line along the right axis. 



(black data points in Fig. [SJa)). Two additional contri- 
butions to the entropy become significant above 0.5 K: 
the phonon contribution, and a magnetic contribution 
arising from thermal population across the spin gap as- 
sociated with the intact dimer triplet states. However, a 
reasonably accurate estimate of the magnetic heat capac- 
ity due to the unpaired Mn ions (ACmag, blue line) can 
be obtained by subtracting the heat capacity of the sto- 
ichiometric parent compound, Ba3Mn208 (red squares). 
The corresponding change in entropy AS* between 50 mK 
and 2 K (right hand axis) is 0.291 J/K mol, substantially 
surpassing S* = 2xi?ln(3/2) = 0.135 J/K mol. 

For comparison, the heat capacity was also measured 
over this temperature range for a crystal with x = 0.980, 
containing approximately the same concentration of un- 
paired Mn moments as the x = 0.020 sample. For 
X = 0.020 the small concentrations of V impurities in- 
troduces unpaired magnetic Mn ions, which interact via 
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the singlet "sea" arising from the majority of unbroken 
dimcrs. In contrast, for x = 0.980, Ba3V208 provides 
an "empty" magnetic background in which the Mn im- 
purities interact. Data are shown in Fig. [5l^b), together 
with an estimate of the magnetic contribution to the heat 
capacity (ACmag, blue line) obtained by subtracting the 
heat capacity of Ba3V2 08 (red squares). ACmag rises 
with decreasing temperature, and appears to begin to 
curve over below approximately 0.5 K. Below 0.2 K a 
clear Schottky anomaly is observed, with the maximum 
centered at 110 mK. The magnitude of the anomaly is 
comparable to that found for x = 0.020, presumably for 
the same reason that a substantial amount of the to- 
tal magnetic entropy is removed at higher temperatures. 
The integrated entropy between 50 mK and 4 K is 0.271 
J/K mol. Measurements to higher temperature arc ham- 
pered by systematic errors introduced by uncertainty in 
measurements of the sample mass for the two crystals 
used in these measurements. However, it is clear that 
the integrated entropy over this window far exceeds that 
which would be estimated if the = ±1 doublet re- 
mains degenerate, S* = 0.135 J/K mol. Hence, for both 
a;=0.020 and 0.980, it appears that interactions between 
unpaired Mn moments leads to progressive removal of 
magnetic entropy over a wide temperature range. And 
in both cases a relatively small number of moments re- 
main to a low enough temperature to contribute to the 
observed Schottky anomaly associated with the single ion 
anisotropy zero field splitting. In neither case is there any 
evidence of a phase transition down to 50 mK. 



IV. DISCUSSION 

In the absence of long range magnetic order, it is rea- 
sonable to consider whether the low temperature state 
of Ba3(Mni_j;V2;)208 is a spin glass, the archetypal 
groundstate for systems with disorder. In canonical 
spin glass systems, the temperature dependence of the 
heat capacity and magnetic entropy scales with moment 
concentration^Si^i. Such systems exhibit a broad fea- 
ture in the heat capacity, the maximum of which occurs 
at a temperature slightly above the freezing tempera- 
ture, with roughly 80% of the total magnetic entropy 
accounted for above Tp. For example, in CuMn, an in- 
crease in the Mn concentration by an order of magnitude 
increases the freezing temperature by greater than a fac- 
tor of 5; such that the magnetic entropy is collected over 
a much wider temperature range. However, the behavior 
in Ba3(Mni_i:Vi:)208 is quite different; in this case the 
entropy is collected over an equivalent temperature range 
independent of vanadium concentration for the range of 
compositions studied (Fig. El^b)). Thus it is unlikely 
that the groundstate is a canonical spin glass. Equally, a 
Griffiths phas o^^'^^ is not anticipated because the parent 
compound itself exhibits a singlet ground state, rather 
than long range magnetic order. 

An alternative groundstate for disordered systems is 



the random singlet phase. This state, which does not 
exhibit long range order, is characterized by formation 
of local singlets^!. Microscopically, as temperature is de- 
creased the pair of most strongly coupled moments forms 
a singlet; then as temperature is further decreased the 
pair of remaining moments with strongest coupling forms 
a singlet, and so on until all the moments arc paired. 
This state was first proposed to account for the low- 
temperature susceptibility and heat capacity of lightly 
doped semiconductors, in particular Si:B^"— , and has 
since been of much theoretical and experimental interest 
for one dimensional disordered systems^"—. The dom- 
inant thermodynamic signature of this phase is the sus- 
ceptibility which diverges as T" with — 1 < a < at low 
temperatures (i.c a sub-Curie law), as observed in both 
Si:P, Cd:S and the one dimensional systems MgTiOB04 
and MnMgBzOg^i^i^ii^. Additionally, it was found 
that the heat capacity of Si:P increases with decreasing 
temperature for the lowest concentrations over an appre- 
ciable range of temperatures reflecting the progressive 
removal of magnetic entropy upon cooling^. 

A scenario in which the ground state of 
Ba3(Mni_a;V2;)208 is a random singlet phase is 
qualitatively consistent with the experimental observa- 
tions. Consideration of the relevant exchange energies 
serves to illustrate whether collective or local behavior 
is expected to dominate. Specifically, a comparison of 
the maximum exchange energy between two impurities 
(JmaxSi ■ Sj) with the average total exchange energy 
resulting from interaction with all of the neighboring 
impurities ZiP^^iJij, where Zi is the number of 
coordinating sites with exchange Jij and Px^i is the 
probability that these sites are occupied) distinguishes 
whether a random singlet phase or collective mag- 
netic order (for example a spin glass) is favored. For 
vanadium concentrations up to the maximum studied, 
X = 0.046, it is straightforward to show that the random 
singlet phase is favored. As an example, for the case 
of x = 0.046 we first consider a free Mn moment for 
which the nearest unpaired Mn is a distance r2 away 
(recall that J2=0.256 meV is the strongest exchange 
interaction after Jq). In comparison, the average sum 
of further interactions acting on the free Mn moment, 
which would be appropriate for a mean field treatment, 
is 3(0.046)Ji + 6(0.046)J3 + 3(0.046)J4 = 0.060 meV. 
The clear energetic advantage of singlet formation 
continues when progressively smaller exchange couplings 
are considered after removal of all J2-bonded pairs, and 
so on (a crude version of the "decimation" procedure 
first introduced to describe random singlet formation3i) . 
This is due to the combined actions of the hierarchy of 
exchange interactions found in Ba3Mn208, the small size 
of the Mn moment {S = 1), and the low concentration 
of magnetic impurities. This crude analysis clearly 
shows that the exchange between a single pair of spins 
dominates over more collective behavior, implying that 
a random singlet state is theoretically favored over a 
spin glass state, at least for the concentration regime 



7 



considered here. 

To further test this view, a model of the random sin- 
glet phase in Ba3(Mni_a;V2;)208 was undertaken in or- 
der to make a quantitative comparison with the mea- 
sured heat capacity data. The model is an approxima- 
tion of the previously described "decimation" process^!, 
with three changes: (1) we do not recalculate the proba- 
bility of finding successively weaker random singlet pairs 
after each round of decimation (this leads to an overesti- 
mate of the probability of smaller pairing distances); (2) 
the exchange interaction between two unpaired moments 
forming a singlet is not renormalized by other exchanges 
between neighboring moments; and (3) we explicitly in- 
clude the effect of single ion anisotropy. Limitations of 
this model are discussed below. The total heat capacity 
is calculated by summing over the probability that two 
moments with exchange J(r) will pair, Px{J{r)), times 
the heat capacity for two moments with exchange J{r): 

CpMT) = 5]F.(J(r))C,(J(r),T) (2) 

J(r) 

Details of the heat capacity calculation for two moments 
with exchange J are given in the Appendix. The ex- 
change between neighboring impurities J{r) was deter- 
mined by the distance between impurities. For neigh- 
boring impurities a distance r < r4 apart, the effective 
exchange was approximated using the values determined 
from INS measurements of the undiluted systemic. For 
impurities at a distance r > r^, the effective exchange was 
approximated by a decaying exponential as expected for 
localized moments: J{r) = J'e"''/'' . The parameters J' 
and r' were determined from fits of the known exchanges 
Ji through J4. The probability Px{J{r)) was deter- 
mined by finding the distance between an impurity and 
the neighboring impurity with which it has the strongest 
exchange, and was calculated similarly to Px{i'), shown 
in Fig. [IJc)^. For r > r4, successive shells of widths 
drn were taken, such that the number of atoms in each 
shell corresponds to a 5% increase in probability (i.e. 
Px{J{rn)) = 0.05 where r„ = r„-i + dr„). 

Both the estimated exchange model and also the prob- 
ability functions used in this model are approximations. 
We have used "bare" exchange values, determined from 
INS measurements, but effective exchange values be- 
tween unpaired Mn moments on broken dimers will be 
renormalized by the presence of the surrounding sin- 
glet states. In particular, geometric frustration naturally 
leads to a reduction in the efi^ective exchange from the 
bare values^. The model further approximates the ex- 
change for r > r4, in which the real effective exchange 
most likely does not vary perfectly exponentially for in- 
termediate r since the superexehange depends sensitively 
on bond angles and bond lengths. Also, Px{J{r)) overes- 
timates the probability of smaller pairing distances. This 
can be understood by considering an impurity i with 
nearest neighboring impurity j. If there is a third im- 
purity k which is closer to impurity j than impurity i, 
then impurity j pairs with impurity k and not impurity 




FIG. 6: (Color online) (Color online) Comparison of the ex- 
perimentally observed heat capacity of Ba3Mn(i_a;V2:)208 
with that of the random singlet model. Panels (a), (c) and (e) 
show C, C /T and the integrated entropy S respectively for the 
random singlet model, as described in the main text. Values 
have been scaled by the vanadium concentration. Panels (b), 
(d) and (e) show the experimental results for the same quan- 
tities. For X = 0.980, the data have been normalized by the 
amount of Mn (i.e. 0.020) rather than the V concentration. 



I, implying that impurity i pairs with a fourth impu- 
rity further away than impurity j, shifting probability 
to longer pairing distances. Nevertheless, as shown be- 
low, this crude model appears to capture the observed 
behavior at a semi-quantitative level, justifying its con- 
sideration in terms of an initial description of the data. 

The calculated heat capacity, the heat capacity scaled 
by temperature and the entropy for the different vana- 
dium concentrations experimentally studied are shown 
in Fig. ini[a), (c) and (e), respectively. For x < 0.046, 
the calculated curves are in reasonable qualitative agree- 
ment with the experimental data, shown for comparison 
in panels (b), (d) and (f) (the experimental data have no 
background subtracted which is negligible in this temper- 
ature range as can be seen in Fig. [S]). For x = 0.009 and 
X = 0.020, the model heat capacity increases as temper- 
ature decreases, reaching a maximum at the single ion 
Schottky peak before decreasing at lower temperatures. 
For X = 0.034 and x = 0.046 the model heat capacity 
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decreases monotonically as temperature decreases, with 
a shoulder at the single ion anisotropy Schottky peak. 

The calculated entropy is also consistent with the ex- 
perimental data (Fig. IHl^e) and (f), respectively) for 
X < 0.046. In spite of the five times difference in 
the vanadium concentration, which leads to a significant 
change in the probability distribution of exchange val- 
ues of the random singlets, the change in entropy from 
50 mK to 500 mK determined from the model rises to 
roughly the same value of ^--^0.38 J/K mol V = 0.42 
-Rln(3) independent of vanadium concentration, in ac- 
cordance with the experimental data. By extending the 
model calculations to both lower and higher tempera- 
tures, the change in entropy for all vanadium concentra- 
tions eventually reaches the maximum available entropy, 
i?ln(3). However, the smaller vanadium concentrations 
recover most of the unaccounted entropy below 50 mK 
while the larger vanadium concentrations recover most 
of the unaccounted entropy above 500 mK. In this sense 
it is somewhat coincidental that the experimentally ac- 
cessible temperature window of 50 - 500 mK leads to an 
almost concentration independent change in entropy. 

Despite the successful description of the evolution of 
the T-dependence of the heat capacity as a function of 
composition, both the calculated heat capacity and en- 
tropy are uniformly smaller than the experimental data, 
reflecting limitations of the simplified exchange model 
used for the calculations. In particular, both the effec- 
tive exchange J(r) and the pairing probability Px{J{r)) 
used in the calculations were approximations as described 
earlier. As a consequence, the model heat capacity de- 
viates slightly from the experimental data. This can be 
seen most clearly by considering Cp/T, which emphasizes 
subtle differences at very low temperatures (Fig. IHKc) and 
(d)). The random singlet model predicts a roll-over be- 
low ~ 80 mK, whereas the experimental data continue 
to rise with decreasing temperature down to our base 
temperature of 50 mK. This difference reflects the in- 
accuracy of the simplified model used to estimate the 
superexchange for the most distant pairs. However, the 
overall agreement is reasonable based on these approx- 
imations, lending considerable weight to the hypothesis 
that Ba3(Mni_a;Va;)208 realizes a random singlet phase 
at low temperature, at least for the range of concentra- 
tions studied in this report. 

It is interesting to compare the heat capacity for x 
= 0.020, corresponding to 2.0% unpaired Mn ions in an 
otherwise "filled" Mn lattice, with that of x = 0.980, 
corresponding to the case of 2.0% Mn impurities in an 
otherwise "empty" lattice. Within the simplified ran- 
dom singlet model developed above, the heat capacity 
for the two cases is identical. However, as can be seen 
in Fig. ini (b), the measured heat capacity for these two 
compositions, although similar, is not identical. This dif- 
ference reflects the renormalization of the bare superex- 
change values for the case of the filled lattice^, and in- 
deed the data for x = 0.980 conform more closely to the 
model (which uses the bare superexchange constants de- 
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FIG. 7: (Color online) (a) Experimental susceptibility of 
Ba3(Mno. 980 Vo. 020)208 for fields both parallel (up triangles) 
and perpendicular (down triangles) to the c axis. Solid black 
line shows Curie contribution to the susceptibility determined 
from a fit of the entire data to Eq. [1] (b) Log-log plot com- 
parison of susceptibility associated with concentration of un- 
paired Mn moments derived from the random singlet model 
(dashed red line) and Curie model (solid black line). Vertical 
dashed line at 1.8 K marks the lowest temperature to which 
experimental measurements reached. For the measured tem- 
perature range, the susceptibility of the random singlet model 
does not deviate substantially from Curie's law. 



rived from INS measurements) than does the data for 
X = 0.020. Nevertheless, both compositions qualitatively 
conform to the expectations of the random singlet model, 
exhibiting the removal of magnetic entropy over a broad 
range of temperatures and a modest feature associated 
with the single ion anisotropy. 

Having introduced the random singlet model for the 
substituted lattice, it is worthwhile to briefly revisit the 
low-temperature susceptibility of Ba3(Mni_a;Va;)208. It 
is straightforward to show that the random singlet model 
results in a sub-Curie power law for the susceptibil- 
ity, as anticipated22i2Si^i^i^i^ (see Appendix). The 
largest exchange interaction between unpaired Mn mo- 
ments (^ J2 ~ 0.256 meV) sets the temperature scale 
below which singlet formation begins to occur. Hence, for 
temperatures above a few Kelvin, the calculated suscep- 
tibility tends towards a simple Curie behavior, as shown 
in Fig. [71[b). For our experiments, susceptibility mea- 
surements were only possible down to 1.8 K, and thus in 
the analysis described earlier we used a simple Curie law 
to fit the data (Eq. [T]). As can be seen from Fig. EKa), the 
measured susceptibility is essentially isotropic over this 
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temperature range, consistent with the small magnitude 
of the single ion anisotropy {D ^ —0.024 meV) and the 
nearly isotropic p-tensor. The observed upturn at low 
temperatures is not perfectly described by the Curie law, 
as anticipated from Fig. [7l[b) , but experimental data ex- 
tending to much lower temperatures are required for a 
quantitative comparison with the random singlet model. 

Finally, we note that although the current data clearly 
indicate the absence of long range order down to 50 mK, 
we cannot rule out an ordering transition at a lower tem- 
perature. However, the relatively small magnetic entropy 
that remains at this temperature, and the progressively 
larger physical separation of unpaired moments due to 
random singlet behavior, make long range order at a 
lower temperature rather unlikely. 



V. CONCLUSIONS 

In summary, the low-temperature properties of the 
site-diluted quantum magnet Ba3(Mni_j;Va;)208 have 
been explored through heat capacity experiments. No 
sharp features associated with a transition into a "order 
by disorder" phase were observed, but rather we find a 
slowly varying magnetic contribution at low temperature 
due to the progressive removal of magnetic entropy over 
an extended temperature range. A Schottky anomaly in- 
duced by the single ion anisotropy is superimposed on 
this behavior below ^--^250 mK. There is very little vari- 
ation in the temperature dependence of the magnetic 
entropy as vanadium concentration is varied, indicating 
that the behavior is not associated with spin freezing. 
Rather, the data suggest that Ba3(Mni_a;Va;)208 has a 
random singlet ground state. 
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The energy spectrum for this system is most easily de- 
termined in the dimer basis. In this case, all the dimer 
states arc diagonal with the first term of the Hamilto- 
nian, JSi • S2. Further, the single ion anisotropy term 
D (('S'l )^ + (Si)^) conserves total and cannot connect 
states of opposite symmetry, implying that only the |00) 
and 1 20) states can have off-diagonal matrix element&i^. 
The matrix elements for all the states are: 
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Appendix A: Derivation of Heat Capacity and 
Susceptibility of a dimer including the effect of 
single ion anisotropy 

Here we numerically calculate the heat capacity and 
susceptibility for a pair of spins with exchange J and 
single ion anisotropy D. This calculation is done for a 
range of antiferromagnetic exchange values relative to the 
easy axis single ion anisotropy, J/\D\. The Hamiltonian 
acting on these two moments is: 



n = jSi ■ S2 + D {{s',f + {s^'f 



•HI 20) 

•HI 00) 

•H|2±2) 
•H|2±l) 
•H|l±l) 
•HI 10) 



3 J 



2D 
T 



|20) 



2^/2D 



100) 



^120) + f 100) 

(3J + 215)12 ±2) 
(3J + i:')|2±l) 
iJ + D)\l±l) 
{J + 2D)\10) 



(A2) 



(Al) 



Diagonalizing the |00) and |22) matrix elements re- 
sults in eigenvectors IS*^) with eigenvalues -B|s±) = 

2D + 3J± ^J{2D + 3jf - WDJ^ . The energy 

spectrum is plotted versus the size of the exchange J/\D\ 
in Fig. m Panel (a) shows the spectrum for large values 
oi J/\D\. In this limit the system approaches the iso- 
lated dimer model, where the singlet has gaps of J and 
3J to the excited triplet and quintuplet states, respec- 
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FIG. 9: (Color online) Heat capacity curves versus normalized 
temperature for specific J/|-D| values plotted on a log scale. 
Mol refers to one mole of spins. 



tively. Panel (b) shows the spectrum for small values of 
J/\D\. In this limit the system approaches the isolated 
moment model, where a dimer has total energy 0, D or 
2D if it is composed of two, one or zero = spins, 
respectively. 

The energy for a pair of moments with exchange J can 
be expressed as a function of temperature: 



Ej{T) 



(A3) 



The summation runs over all the different basis states 
(|7,) = |5±),|2±2),|2±l),|l±l),|10)),i;|^,) is the en- 
ergy of a given state and /? = l/ksT. The heat capacity 
of a pair of moments with exchange J is easily calculated 
numerically from the energy: 



Cp{J,T) 



d{E.j{T)) 
dT 



(A4) 



Several different heat capacity curves for a range oi J /\D\ 
values are plotted on a log scale in Fig. [S] For J ~ Q 
(solid black line) , the system reduces to the isolated mo- 
ment case and the heat capacity has a peak induced by 
the two gaps of the single ion anisotropy (a gap of D, 
2D between the \S~), |10), |2 ± 2) ground states and the 
|2 ± 1), |1 ± 1) first excited states and \S^) second ex- 
cited state, respectively). For the largest J values of 



J = SlDljlOlDj and 30|Z)| there is a single peak cen- 
tered at roughly 0.4 J arising from both the singlet-triplet 
and singlet-quintuplet gaps. For J ~ \D\ there is a sin- 
gle peak with a shoulder at lower temperatures. For 
J = 0.3|Z)| there are split peaks, one centered roughly at 
the same position as the single ion anisotropy peak and 
a second at lower temperature. Finally for the small- 
est non-zero J values of J = 0.1|D| and 0.03IZ?! there 
are three peaks stemming from three gaps from the \S~) 
groundstate: the gap to the first excited state, |10); the 
gap to the second excited states |2 ± 2); and finally the 
single ion anisotropy gap to the third excited states, 
|1 ± 1), |2 ± 1) and \S+). The model shown in Fig. El 
used a superposition of heat capacity curves taken from 
exchange values mostly within this range. The different 
model curves used differently weighted superpositions, 
such that the higher (lower) vanadium concentrations are 
weighted more towards larger (smaller) J values based on 
the probability distributions Px{J{r)). 

The susceptibility is calculated in the low field limit by 
similar means. The calculation was done with H\\c (i.e. 
along the dimer direction) for simplicity since most of 
the dimer states described above are diagonal and there 
is no appreciable anisotropy based on field direction for 
the temperature range studied (Fig. [D^a)). The low field 
susceptibility for a pair of moments with exchange J as 
a function of temperature is: 



M{J,T) 1 9l^BT.,Tn'^^^f 



H 



H 



(A5) 



17.) 



is the z component of the angular momentum of 
energy eigenvector [7^), appropriate for fields along the c 
axis. The total susceptibility is calculated by summing 
over the probability that two moments with exchange 
J(r) will pair, Px{r), times the susceptibility of two mo- 
ments with exchange J{r): 



Xtot 



{T) = Y,Px{J{T))x{J{r),T) 



(A6) 



J{r) 



The susceptibility calculated from this random singlet 
model for an a; 0.020 system (dashed red line) is plotted 
in a log-log scale in Fig. [71[b) . 
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